Fidelity analysis of topological quantum phase transitions 
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We apply the fidelity metric approach to analyze two recently introduced models that exhibit a quantum 
phase transition to a topologically ordered phase. These quantum models have a known connection to classical 
statistical mechanical models; we exploit this mapping to obtain the scaling of the fidelity metric tensor near 
criticality. The topological phase transitions manifest themselves in divergences of the fidelity metric across 
the phase boundaries. These results provide evidence that the fidelity approach is a valuable tool to investigate 
novel phases lacking a clear characterization in terms of local order parameters. 
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Introduction. — This is an exciting period for condensed 
matter physics, when novel phases of matter that defy tradi- 
tional understanding are being observed and predicted. Ex- 
amples include topological phases [ 1 ] which cannot be de- 
scribed by Landau-Ginzburg- Wilson paradigm Absence 
of local order parameters and symmetry breaking mechanisms 
are among the most remarkable features of these systems. 
These novel phases arise, for example, in collective phenom- 
ena exhibited in strongly correlated systems of two dimen- 
sional electrons at very low temperature, like in the fractional 
quantum Hall effect yl 01 . In such systems, the motion of 
electrons is highly constrained, and the fluctuations are en- 
tirely quantum in nature. In this situation Landau's theory, 
which is essentially a theory of classical order, can fail. 

It is compelling to find new ways to analyze such phases. 
Using tools from quantum information, it has been possible 
to characterize topological order using the concept of topo- 
logical entropy j5[]-0]. Here, we call for a new information- 
theoretic tool for studying quantum phase transitions (QPTs) 
^ to topological phases. The new notion is the fidelity of 
ground states, whose role in the study of QPTs has been de- 
veloped in Irjjl- lO . The basic idea is that near a quantum 
critical point there is a drastic enhancement in the degree of 
distinguishability between two ground states, corresponding 
to slightly different values of the parameter space that defines 
the hamiltonian. This distinguishability can be quantified by 
the fidelity, which for pure states reduces to the amplitude of 
inner product or overlap. This approach is suitable for detect- 
ing QPTs and analyzing topological phases, since the method 
does not rely on constructing an order parameter, nor on the 
symmetries of the system. The overlap of two nearby ground 
states is a global quantity of the system that does not depend 
on local features like the existence of a local order parameter. 
Therefore, it should contain all the information that describes 
topological order. The capability of fidelity to spot a topo- 
logical QPT has been shown in 112211 by numerical analysis. 
Since topological order is a property of the ground state wave- 
function alone, knowledge of the ground state of the system is 
sufficient in order to carry out this analysis. 



In this work we analyze two models that exhibit topologi- 
cal order: hamiltonians with a stochastic matrix form 113411 and 
the quantum eight-vertex model 13511 . Both present a transi- 
tion from a non-topologically ordered phase to a topological 
phase. Moreover, they exhibit a close connection to classical 
statistical models. Indeed, the fidelity and its second deriva- 
tive, the so-called fidelity metric, are related to the partition 
function ll36ll and correlation functions of the corresponding 
classical model, respectively. 

Topological QPT in stochastic matrix form hamiltonians. — 
In this section we apply the fidelity approach to analyze the 
quantum phase transition to a topologically ordered phase for 
a particular case of hamiltonians that exhibit a stochastic ma- 
trix form decomposition lf3~4h - In the authors showed their 
model had a transition from a magnetically ordered state to a 
topological ordered phase, with the topological entropy hav- 
ing a jump to a nonzero value at the transition. We now ap- 
ply the fidelity approach to this topological QPT and find the 
scaling of the fidelity metric near criticality. Let us start by 
briefly reviewing the model. Given a square lattice with pe- 
riodic boundary conditions and spins 1/2 on the bonds, con- 
sider the following hamiltonian: 
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where A s — Y\ ies of and B p = Y\ i£ of are the star and pi 



quette operators of the Kitaev model 113 811 . The hamiltonian 
TJ with Ao = is said to be written in stochastic matrix form 
33]. 

The ground state of this hamiltonian can be computed ex- 
actly and is given by 1137m 1 : 
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Here, G is the abelian group generated by all the star oper- 
ators A s , |0) is the completely polarized state corresponding 
to all the spins in the +1 eigenstate of a z and crf(g) is the 
z-component of the spin at site i in the state g\0). 

Let us try to understand the phases of this model. When 
(3 = 0, we have the pure Kitaev toric code. Its ground state 
is a closed string condensed phase. An x(z)— string is a col- 
lection of spins that are flipped in the a z (a x ) basis. The term 
with the plaquette operator says that only closed strings are 
allowed. The term containing the star operators A s makes 
instead closed strings of flipped spins to be created and fluc- 
tuate. This phase is topologically ordered, as is shown by a 
non- vanishing topological entropy |[5|]-||3]. We can regard the 
/3-dependent term as a kind of tension for the z— strings. As 
we increase (3, larger loops are less favored. Indeed one can 
see that for small (3 the model is the toric code in an exter- 
nal magnetic field. For larger (3 the phase is not topologically 
ordered, as one can infer from the vanishing of the topologi- 
cal entropy [37]. This is why one can use topological entropy 
as an order parameter 1 22 , 37[ S . One expects that for a 
particular value of (3 the system undergoes a QPT from the 
the topologically ordered phase to a "magnetically" ordered 
phase. The authors in [37] proved that this model has a second 
order phase transition at (3 C = (1/2) ]n(y/2 +1). For (3 < (3 C 
the system has a topologically ordered phase, with Stopo = 2, 
and for (3 > (3 C the topological entropy vanishes: Stopo = 0. 
It is very important to notice that despite being not topolog- 
ically ordered, such phase is not a Landau-Ginzburg phase. 
There is no local order parameter to characterize it 13711 . 

We now analyze this transition using the fidelity between 
two ground states \(3) and \(3 + 8(3) corresponding to slightly 
different values of the relevant parameter (3. Therefore, we 
consider this quantity: 
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Expanding to second order in 5(3 i.e., F w 1 — gpp5(3 2 , 
we obtain the following^ffe/ify metric II 1 3fl gpp: 
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Much in the spirit of the fidelity approach, near the quan- 
tum phase transition there is an enhancement in the distin- 
guishability between the ground states \(3) and \(3 + 5(3), re- 
sulting in a superextensive scaling of the singular behavior 
of this metric at the critical point. Indeed, this singular be- 
havior can be captured by mapping this quantum model to a 
classical statistical model, in the following way. Any group 



element g G G is the product of the star operators in some 
set S(g). Thus, g\Q) is completely specified by the same 
set, modulo the product of all such operators, equal to the 
identity for periodic boundary conditions (i.e., for a torus of 
genus one). Then, for every two configurations specified by 
{g E G} there will correspond one configuration {9} of a 
classical Ising model with degrees of freedom 9 S on the sites, 
such that 6 S = — 1(+1) when the corresponding star op- 
erator A s is (is not) acting on the site s. Since a spin erf 
can be flipped only by its two neighboring 9 spins, we have 
that &i — 9 S 9' S , with s and s' the end points of the bond 
i. In that case, defining _Ei s i ng = JJ2< S s'> ^s0' s , we ob- 
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2Z{(3), where we took (3 = J/T for the Ising model. Using 
this equality, we can write (O as: 
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where C v is the specific heat of the 2D Ising model. It is well 
known that C v has a logarithmic divergence at criticality 114 111 . 
Hence mapping to the classical Ising model reveals that the 
fidelity metric has a logarithmic divergence 
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at/3 c = iln(V2+l). 



In 113711 the authors remark that indeed the phase transi- 
tion to the topologically ordered phase could be detected by 
the local magnetization m((3) = ^J2i(&i) = jj E isin g ((3), 
with its first derivative equal to the specific heat, i.e., ^ = 
— jff3 2 Ci s i ng ((3), where N is the number of sites. We see 
that the fidelity metric captures very naturally this divergence 
since it is equivalent to the specific heat, which diverges at the 
critical point. 

Topological QPT in the quantum eight-vertex model. — We 
now turn to analyze another model that exhibits a transition to 
a topological phase, and in which the mapping to a classical 
statistical model can be performed to analyze the scaling of 
the fidelity metric near the critical point. This model is the 
so-called quantum eight-vertex model, defined and studied in 
references 13511 and 14011 . We proceed to review this model 
very briefly. 

The classical eight-vertex model 14 111 consists of arrows 
placed along the bonds of a square lattice. The arrows can 
point in either direction along each of the bonds, subject to the 
constraint that an even number of arrows go into (and out of) 
each site. There are eight distinct configurations for the arrows 
around each site satisfying this constraint. Each vertex con- 
figuration is assigned an energy e^. Furthermore, by imposing 
toroidal boundary conditions, symmetry under rotations and 
inversions of all spins (i.e., zero external electric field), one 
finds that there are only two independent Boltzmann weights, 
usually denoted by c and d, with c = e~ tc -/ T and d = e~ Cd / T . 
The partition function of this model has the form Z(c, d) = 
^2 C c ™c(C)^n d (c)^ w j m Uc (C) and rid{C) the number of c and 
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d type vertices for the configuration C. The total energy for a 
configuration C is given by E — n c (C)e c + nd(C)ed- 

This classical model can be exactly solved in the thermo- 
dynamic limit by computing the free energy density using the 
highest eigenvalue of the transfer matrix 14111 . It exhibits or- 
dered phases for d > c + 2 and d < c — 2, where the Z 2 
symmetry of flipping all the arrows is spontaneously broken, 
while the system is disordered for |c — d\ < 2. For d < c — 2, 
there is a proliferation of c vertices, while for d > c + 2 the d 
vertices dominate. These phases are called "antiferroelectric". 

One unusual feature of the classical model is that the crit- 
ical exponents change continuously along the critical lines 
d = c + 2 and d = c — 2, with the free energy density having 
a singular behavior near d = c — 2 of the form 
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with fi = 2tan~ 1 \fcd. When tt/ [i — to, with m an integer, 
this expression is changed by an additional logarithmic diver- 
gence: / s i n g ~ | | c ;_ c |_2| 7r /^ ln(|d— c|— 2). The model is also 
critical along the lines c = 0, d < 2 and d — 0, c < 2, since it 
reduces to the disordered phase of the six-vertex model, which 
has an infinite correlation length. The points c = 0, d = 2 and 
d = 2, c = are BKT critical points. There, the exponent 
7r//i diverges. 



The quantum eight-vertex model [35] is defined such that 
its Hilbert space basis {\C)} is given by the configuration 
space of the classical eight-vertex model, with each state real 
and orthonormal to each other. The hamiltonian of this model 
is of the form H = ^\ Qi, with Qi positive operators, chosen 
such that H annihilates the following state: 
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operators for the c and d type vertices, for the configuration C 
IHoh . The authors in Q and H noted that since the normal- 
ization factor above is the partition function for the classical 
two-dimensional eight- vertex model with weights c 2 and d 2 , 
then the ground-state phase diagram for the quantum model is 
identical to the classical one, but given in terms of c 2 and d 2 . 
The quantum model exhibits a topologically ordered phase in 
the region of the phase diagram that corresponds to the dis- 
ordered phase of the classical model. Indeed, the topological 
entropy in the quantum model is given by Stopo = — m (2) in 
the topological phase \d 2 — c 2 \ < 2, while it is zero elsewhere. 
In particular, for c 2 = d 2 = 1 one recovers the ground state 
of the Kitaev model H]. 

Let us now pursue a fidelity analysis of this quantum phase 
transition. Again the mapping to the classical model proves 
useful. As we will see, the fidelity metric is equal to the fluc- 
tuations in the number of c and d type vertices of the classical 
model. This will provide us with the scaling of the metric near 
the phase transition. 



Consider then the fidelity between two ground states for 
slightly different values of the parameters c 2 and d 2 and ex- 
pand it to second order in c 2 andeP: F = (gs(c 2 , d 2 )\gs(c 2 + 
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where the averages are now taken with respect to the classical 
eight- vertex model. 

Using this equivalence with the classical model, we can 
get the scaling of those metric elements near criticality 
by using the expression for the free energy density / = 
— TliniAT^oo iV _1 In Z(c 2 , d 2 ) as a generating function for 
correlations, by differentiating with respect to the energies e c 
and Ed- We obtain then that the dominant scaling near critical- 
ity of the metric elements is: 
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Then, we have an algebraic divergence of the fidelity metric 
only for n/fi — 2 < 0, and a logarithmic divergence for 7r//j — 
2 = 0. Using the fact that near criticality /i = 2 tan -1 v c 2 d 2 , 
those two conditions can be written as 1 < c 2 d 2 and c 2 d 2 = 1, 
respectively. Contrary to the case analyzed before, the metric 
now diverges as a power law instead of logarithmically, but 
only for a certain region of the phase diagram. 

Some remarks are now due. The eight-vertex model can be 
shown to be equivalent to two classical square Ising lattices, 
coupled with a quartic spin term [41], with both models at 
the same temperature. It is interesting to note that the curve 
c 2 d 2 — 1 corresponds to the line along which the coupling 
between the four spins disappears, and separates the region 
where this coupling is ferromagnetic and anti-ferromagnetic. 
The region 1 < c 2 d 2 corresponds to this last case. 

Furthermore, given the equivalence of the eight vertex 
model to a pair of decoupled Ising models for c 2 d 2 = 1, one 
can contemplate the effects of disorder onto the previous scal- 
ing results by considering the two dimensional Ising lattice 
model with random bonds. For c 2 d 2 = 1 the metric elements 
(t9T>-(fTTT> are now proportional to the specific heat C v of the 
Ising 2D model, which can scale as log(logi) at criticality, 
with t = (T — T c )/T c the reduced temperature JH, iH . 
We conclude that the fidelity metric elements can present this 
novel doubly logarithmic divergence for the disordered quan- 
tum eight vertex model as well. 

Finally, there are many models which are either a special 
case of the eight-vertex model or equivalent to it, such as the 
quantum XYZ chain model, ice model, F model, etc 14111 . 
Therefore, by analyzing the fidelity metric in the quantum 
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eight-vertex model we can infer its behavior in many other 
models and predict its divergence only for regions in the 
phase diagram that correspond to the condition 1 < c 2 <P. 

Conclusions and outlook. — In this paper we have per- 
formed a fidelity analysis of quantum phase transitions to 
topologically ordered phases. We have considered two dif- 
ferent systems that can be naturally mapped onto classical 
statistical mechanical models. The mapping reveals that the 
fidelity metric corresponds to derivatives of the free energy 
with respect to some parameters of the model, giving rise to 
correlations in the classical system. 

We discovered a logarithmic divergence in the fidelity met- 
ric for the stochastic matrix form Hamiltonian model near the 
transition to the topologically ordered state. This may be re- 
lated to the fact that in this model there is no local oder pa- 
rameter nor symmetry breaking, and only topological order is 
involved. In perspective, this is an aspect that deserves more 
investigation. On the other hand, the quantum eight vertex 
model still has a power law divergence at the transition to the 
topological phase, but exhibits this singularity for a restricted 
region of the phase diagram only. 

A satisfactory understanding of the relation bewteen fidelity 
metric singularities and the nature of different topological as 
well as standard orders involved in the transitions is a primary 
goal for future investigations. 
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Note added. — Upon the completion of this work, we be- 
came aware of related results for topological phases in the Ki- 
taev honeycomb model in H44I1 and 14511 . 
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